Preprint typeset in JHEP style. - HYPER VERSION 

SLAC-PUB-8976 
UCLA/Ol/TEP/16 
September, 2001 



Two-Loop Amplitudes for Gluon Fusion into 
Two Photons 



Z. Bern,* A. De Freitas* 

Department of Physics and Astronomy 
UCLA, Los Angeles, CA 90095-154? 

L. Dixon^ 

Stanford Linear Accelerator Center 
Stanford University 
Stanford, CA 94309 



Abstract: We present the two-loop matrix elements for the scattering of two gluons 
into two photons in QCD. These matrix elements will enter into improved estimates of 
the QCD background to Higgs production at the LHC when the Higgs decays into two 
photons. The photon mode is especially important if < 140 GeV. We also give the 
amplitudes for the crossed process, (77 — > g'j. 



Keywords: |two-loop, QCD, quantum chromodynamics, Higgs, photons| . 



Submitted to JHEP 



*Research supported by the US Department of Energy under grant DE-FG03-91ER40662. 
'''Research supported by the US Department of Energy under grant DE-AC03-76SF00515. 



1. Introduction 



One of the most pressing problems in particle physics today is to determine the nature 
of electroweak symmetry breaking. Experiments in the coming decade at the Fermilab 
Tevatron and the CERN Large Hadron Collider (LHC) will address this problem, in 
particular by searching for one or more Higgs bosons. There are good reasons to think 
that at least one such particle will be fairly light. In the Standard Model the Higgs 
boson is constrained to be light by its influence on precision electroweak measurements, 
ttih ~ 205-230 GeV at 95% CL [|I|. In the Minimal Supersymmetric Standard Model, 
the lightest Higgs boson is predicted to have a mass below about 135 GeV [0]; over 
much of the parameter space it has properties reasonably similar to the Standard Model 
Higgs boson. Finally, there are hints from LEP at 2.9a significance of a Higgs boson 
at ruH = 115 GeV [|. 

To aid in the search for the Higgs boson, it is useful to have a detailed understand- 
ing of the Standard Model backgrounds. At the LHC the most important mode for 
discovering a Higgs boson with niH < 140 GeV is via its decay into two photons ^. 
The irreducible two-photon background due to QCD, pp — > jjX, will be determined 
experimentally at the LHC, so the Higgs search can proceed even with imprecise the- 
oretical knowledge of it. Nevertheless, it is still of interest to have robust theoretical 
predictions prior to the experiments to help optimize Higgs search strategies. In this 
paper we provide the two-loop matrix elements for gg —>■ 77 with massless quarks in 
the loop, which form a central ingredient in an improved prediction of the irreducible 
background. 

The process pp — *■ 77X proceeds at lowest order via the partonic subprocess 
qq 77, which is independent of a^. The next-to-leading-order (NLO) corrections 
to this subprocess have been incorporated into a number of Monte Carlo programs 0. 
However, at the LHC, due to the high gluon luminosity, formally higher order correc- 
tions involving gluon initial states are sizable. The contribution from the lowest order 
gg — > 77 amplitudes (from one-loop box diagrams), though of order a^, is very compa- 
rable in size to the qq 77 contribution @, 0, ^ . Thus a calculation of the gg — > 77 
subprocess at its next-to-leading-order (which is formally N^LO as far as the whole pro- 
cess pp — > 77X is concerned), should lead to a significant reduction in the uncertainty 
on the total cross section. A calculation of the background, incorporating the two-loop 
matrix elements presented in this paper, will appear in a forthcoming publication 
In the range of di-photon invariant masses relevant for the Higgs search, 90-150 GeV, 
quark masses may be neglected. The u, d, s, c and b quark masses are all much less 
than the scale of the process, while the top quark contribution is negligible until the 
invariant mass approaches 2mt ~ 350 GeV. 
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The subprocess (77 g^y does not possess quite the same phenomenological signif- 
icance as gg 77. However, it does receive a large power-enhanced correction at two 
loops in the forward scattering limit, due the exchange of two gluons, which represents 
the beginning of the Reggeization of this amplitude. The amplitudes for gg —>■ 77 and 
91 ~^ 91 cire of course related by crossing symmetry, but writing the full two-loop am- 
plitude in a crossing-symmetric form would be more cumbersome, so we shall present 
each case separately. 

Calculating two-loop four-point amplitudes involving more than a single kinematic 
variable is a relatively new art. The first calculations of this type in gauge theory were 
for the special cases of gluon-gluon scattering with maximal supersymmetry ||^ and 
maximal helicity violation in QCD |jlO|. More recently, more general calculations of 
interferences of two-loop amplitudes with tree amplitudes in QED and in QCD have 
appeared. In QED the interferences for e~^e~ fJ'^f^~ and e+e~ —>■ e+e~ have 
been computed. In a tour de force calculation, Anastasiou, Glover, Oleari, and Tejeda- 
Yeomans have provided the interferences for all QCD 2 2 parton processes |T^, |T^ . 

The two-loop Feynman diagrams for the gg 77 matrix elements presented in 



this paper are similar to those required for gluon-gluon scattering [0, gg ^ gg, except 



that many of the non-Abelian diagrams are not present. In the gg 77 case, the tree 
amplitudes vanish and the one-loop amplitudes give the leading order contributions. 
Thus the next-to-leading order contributions to this process require a different inter- 
ference, of two-loop amphtudes with one-loop amplitudes. Instead of evaluating this 
interference directly, we have computed the two-loop gg 77 amplitudes in a helic- 
ity basis. We employed a unitarity- or cut-based technique Jl^ ||, to generate the 
required loop momentum integrals. These integrals were then evaluated using recently 
developed techniques. 

Two important technical breakthroughs which have provided the required integrals 
for general 2 — >■ 2 scattering in the massless case are the calculations of the dimen- 
sionally regularized scalar double box integrals with planar [jl5| and non-planar 
topologies and all external legs massless, and the development of reduction algorithms 
for the same types of integrals with loop momenta in the numerator (tensor inte- 
grals) |]T^ |T8|, |T^. Related integrals which also arise in the reduction procedure have 
been computed in refs. Taken together, these results are sufficient to compute all 
loop integrals required for 2 ^ 2 massless scattering amplitudes at two loops. 

A helpful development for performing explicit two-loop calculations is a general 
formula due to Catani for the infrared divergence appearing in any two-loop QCD 
amphtude |2ll]. By appropriately adjusting group theory factors, it is straightforward 
to apply Catani's QCD formula to amphtudes including external photons. Because of 
the vanishing of the tree amplitudes, the infrared divergences for gg — > 77 and (77 g'^ 
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are much tamer than those of a typical QCD process at two loops. Indeed, Catani's 
formula collapses to a form previously derived for one loop amplitudes p2|, We 



present the amplitudes in terms of Catani's formula for the infrared divergences, plus 
finite remainders for each independent helicity configuration. 

In the physical next-to-leading-order correction to the gg 77 subprocess of 
PP — > 77X, the infrared divergences from the two-loop amplitude cancel against those 
arising from phase-space integration of the (square of) the one- loop amplitudes for 
gg — i> 775f [p4|, p5| , after factorizing the initial-state collinear singularities. Both these 
"virtual" and "real" divergences are effectively the same as encountered at one loop. 



so standard one- loop formalisms ||26| can be employed to obtain the cross section 0]. 

This paper is organized as follows. In section ^ we review the infrared structure 
of the two-loop gg 77 and (77 (?7 amplitudes. The corresponding one-loop 
amplitudes, which appear in the formula for the infrared divergences of the two-loop 
amplitudes, are presented in section ^. In sections | and |^ we present analytic results 
for the finite remainders of the two-loop gg 77 and (77 g'y helicity amplitudes, 
respectively. In section ^ we discuss the checks we performed on our results. Finally, 
in section |^ we present our conclusions. 



2. Review of infrared structure 



In presenting the amplitudes it is convenient to separate the infrared divergent parts 
from the finite parts. Dimensionally regulated two-loop amplitudes for four massless 
particles in D = 4 — 2e dimensions generically contain poles in e up to Catani 
has presented a general formula for the structure of infrared divergences of any QCD 
amplitude |2lJ. With minor adjustments, Catani's formula is also valid for the case of 
mixed amplitudes of QED and QCD. We shall therefore adopt his notation in presenting 
our results. 

We work with ultraviolet renormalized amplitudes, and employ the MS running 
coupling for QCD, ^^(/i^). Since the tree amplitudes for the process under consideration 
vanish, for purposes of renormalization we only need the one-loop relation between the 
bare coupling a" and renormalized coupling as{fi'^), 



(2.1) 



where = exp[e(ln47r + V'(l))] and 7 = -V^(l) = 0.5772 . 
coefficient appearing in the QCD beta function is 

IICa - 4:TRNf 



12n 



is Euler's constant. The 



(2.2) 
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where Nf is the number of hght (massless) quarks, Ca = N for SU{N) and Tr = 1/2 
for fundamental representation Dirac fermions. 

The renormahzed QCD corrections to the gg —>■ 77 amphtude discussed in this 
paper may be expanded as 



M 



aa-*ii 



{p}) = Ana 



{P}) (2.3) 
^^VaiSU-,(a'; {?}) + 0(aUl^')) 



where {p}) is the Lth loop contribution. Since there is no direct coupling 

between gluons and photons, the expansion starts at one loop. 

The QCD MS counterterm which is subtracted from the bare two-loop amplitude 



to obtain the renormahzed -M^^g- 



■77 



(/x^; {p}) is (c/. eq. (O 



C.T. 



IIN - 2Nf 1 
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(2.4) 



The relative simplicity of the ultraviolet subtraction term is due to the vanishing of the 
tree-level amplitudes. 

There is also a QED divergence associated with the dia- 
grams in Figure |I[ In dimensional regularization with mass- 
less fermions, these diagrams would vanish by virtue of con- 
taining a scale free integral. This vanishing represents a can- 
cellation between infrared and ultraviolet divergences. In 
QCD with MS renormalization such diagrams (with the bub- 
ble on a gluon leg) are thus taken to vanish. However, if 
one renormalizes QED in the conventional on-shell scheme, to 
avoid the infrared divergences one should retain the fermion 
masses in the external bubbles. Now the bubble integral is 
nonzero and ultraviolet divergent. But this divergence, and 
indeed the entire integral, is exactly canceled by the on-shell- 

scheme counterterm, precisely because the external leg is a real, on-shell photon. In 
the on-shell scheme, the coupling constant should of course be set to a = a{0) = 
1/137.036. . .. This value should then be used for all the QED couplings associated 
with real, external photons. 

The infrared divergences of a renormahzed two-loop amplitude, for the case where 



%4 



Figure 1: These dia- 
grams give QED diver- 
gences. They are exactly 
canceled by conventional 
on-shell renormalization. 



the tree amplitude vanishes, are [Bl 



IM'H/^'; W))Ka = {p}) \M^^\^'; {p})U, + IM^^^^'^'i^'; {p})U. , (2.5) 
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where IM^^^fi"^; {p})) Rs. is a color space vector representing the renormahzed L loop 
amplitude. The subscript rs stands for the choice of renormalization scheme, and fi is 
the renormalization scale. These color space vectors give the amplitudes via, 



M„(r\...,n'^") = {ai,...,an\MniPi,.-.,Pn)) 



(2.6) 



where the a, are color indices. The divergences of are encoded in the color operators 

/«(6,/.^M). 

In pure QCD, the operator I^^\e,fi'^] {p}) is given by 



l''\e,^';{p}) = l^^,j:T.T^-T, 



2r(l-e) 



i=l jj^i 



1 : 

e2 ^ T? e 



2pi ■ Pj 



(2.7) 



where \ij = +1 if z and j are both incoming or outgoing partons and = otherwise. 
The color charge Tj = {T"'} is a vector with respect to the generator label a, and an 
SU{N) matrix with respect to the color indices of the outgoing parton i. For external 
gluons T^^fe = ^/™^ so T^^ =Ca = N, and 
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The two processes we consider here are 

5((-pi,-Ai) + 5'(-P2,-A2) 
gi-pi,-\i) + 7(-P2,-A2) 



liP3,h) + 7(P4,A4), 
^(^3,^3) + 7(P4,A4), 



(2.8) 

(2.9) 
(2.10) 



using an "all-outgoing" convention for the momentum {pi) and helicity (Aj) labeling. 
The Mandelstam variables are s = {pi + P2Y, t = {pi + PiY, and u = {pi + psY- To 
apply eq. ( |2.7| ) to these processes, one may convert two of the gluon legs to photons by 
setting ■ Tj — > 0, Tj ■ ^ to obtain the simplified formula 



-N 



r(l-e) 



6 3 / e 



(2.11) 



The color structure of the two-gluon, two-photon amplitude can only be proportional 
to 5"^ Hence /^^^_^^^ is necessarily proportional to the identity operator, and the 
color space language is actually unnecessary. In the s-channel where s > one should 
use the usual analytic continuation 



-s -le) 



-e^iire ©(s) 



(2.12) 



For the crossed process g'j gj, the same formulae apply with the obvious modifica- 
tion. 



/liL„.(e,/i';M) 



r(l-e) 



11 liV/M 



^ ^ ' 6 3 N ) e 



-u 



(2.13) 
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The renormalized two-loop amplitudes may thus be separated into divergent parts 
(which also absorb some finite terms) and a finite remainder, 

•^£^-77 = 4^77(^' /^'; {p})Mf^-.,, + Mf^%^ , (2.14) 
^^7^7 = /^'; {P})M^^^,, + M^X, . (2.15) 

The finite remainders Adf^^^^ and A^^^-*5g7 will be presented in sections ^ and |^. 
Since /*^^^ contains 1/e^ poles, the above decomposition requires the expansion of M.^^^ 
through O(e^), which will be given in section]^. 

We present the two-loop amplitudes for definite external helicities, using 't Hooft- 
Veltman (HV) dimensional regularization P^. At one loop, scheme conversions be- 



tween the various fiavors of dimensional regularization have been extensively discussed pS 
Following a similar strategy, one may convert our results in the 't Hooft-Veltman scheme 
to the conventional dimensional regularization (CDR) scheme. In the CDR scheme one 
would compute the interference of say A^^^)_^^^ with A^^y_^^^, summed over all colors 
and helicities. This CDR interference may be obtained from our HV result by using 
eq. ( p.l4|) for Mf^^^^ and replacing the HV interference Mf^^^^ x M^g^Z'yy by the 
same quantity evaluated in the CDR scheme, while leaving J*-^^ and the finite remain- 
der Aif^^^^ of section || unaltered. The same considerations apply for — >• as 
well. When constructing the cross section for pp — >■ 77X, the difference between the 
HV and the CDR scheme amplitudes is cancelled by a similar scheme dependence in 
the real emission terms, and at the end only the finite remainder contributes from the 
two-loop amplitudes [§]. 



3. One- loop amplitudes 



As noted above, the decompositions ( p.l4| ) and ( p.l5| ) require the series expansions of 
the one-loop gg — > 77 and (77 — *• 5^7 amplitudes through C(e^), A simple way to obtain 
these expansions is to express them in terms of four-gluon amplitudes which are known 
to all orders in e p9|, in terms of integral functions whose series expansions are known 
to the appropriate order fT^, 0. The fermion loop contributions to gg — > gg satisfy 
the color decomposition 



A^(i)/ ( 

99^99^ 



1,2,3,4:) = Nf J2 tr[T"-wr''-(2)TM3)TM4)] 
X AW(a(l),a(2),a(3),a(4)) 



(3.1) 



where is the number of quark fiavors and the sum of permutations a runs over 
the six non-cyclic permutations of external legs. We have adjusted the normalization 
factors in the amplitudes to match the conventions of ref. . 
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One may obtain amplitudes with photons from eq. (|3.1|) simply by replacing the 
appropriate color matrices with the identity matrix and altering the normalizations. 
In particular, the amplitudes where any two legs are gluons and the other two legs are 
photons are all the same, 



Nf 



4 = 1 



where Qi are the electric charges (2/3 for up-type quarks and —1/3 for down-type 
quarks) and 

M« = Y: A^'\cr{l), a(2), a(3), a(4)) . (3.3) 

The permutations a run over the same six orderings of external legs as in eq. ( p.l| ) , while 
the overall factor of 2 arises from our non-standard normalization of color matrices, 
requiring an extra factor of for each gluon converted to a photon. In the permutation 
sum the ultraviolet and infrared divergences appearing in the four-gluon amplitudes 
cancel, leaving a finite expression. 

The amplitudes are evaluated in the spinor helicity formalism . It is convenient 
to extract overall spinor phases from each helicity amplitude, 

/l(l)|'lAl r,A2 oAa aX4\ _ C 

/i [l ,Z ,6 ,4 J — OA1A2A3A4 ^AiA2A3A4 ' 

M«(l'^2'^3^^4^^) = S,,x,^,x,Mi^X^^^^^, (3.4) 
where the Aj signify the helicities of each leg and 

S,,,, = ,Ji4M, (3.5) 

++++ (12) (3 4)' ^ ^ 

^ . (12) (14) [2 4] ^ , (2 3) (2 4) [3 4] ^ , (3 2) (3 4) [2 4] 
+++ (3 4) (2 3) (2 4) ' (14) (31) (3 4) ' (1 4) (2 1) (2 4) ' 

, (12) [3 4] _, (13) [24] _, (23) [14] 
[12](3 4)' [13] (2 4)' [2 3] (14)- 

Our conventions here are that all external legs are outgoing; for incoming legs one should 
reverse the helicities of those legs. The spinor inner products [jS^, ^ are {ij) = 



and [ij] = where are massless Weyl spinors of momentum fcj, labeled with 

the sign of the helicity. They are anti-symmetric, with norm \ (ij) \ = \ [ij] \ = ^/SiJ, 
where Sij = 2ki ■ kj. From these definitions it follows that the 5'aiA2A3A4 are phases. 

The fermion loop contributions to the one-loop four-gluon helicity amplitudes, in 
a form valid to all orders in e (with four-dimensional external momenta pi) are 

= e(l-e) Box(^)(s,t), 



7 



A 



(1) 



su 
t — u 



tu 
s — u 



St 



, eBub(^)(t) + — eBox(6)(s,t) 
2u 



+ e(l-e)Box(^)(s,t) 



A 



A 



(1) 

++- 

(1) 



A 



(1) 

- + + + 5 



-iseBox('^)(s,t) - ^Bub(^)(t) 

+ ^Bub^^) (t) - - e Bub(^) (t) + e(l - e) Box^^) {s, t) 



A 



(1) 



^tBub(^)(s) + isBub(^)(i) - ^^^(1 - e) Box(^)(s,t) 



u 

St 

u u u 

_ f^Bub^^H^) - -Bub(^)(t) - -Bub^^H^) + -eBub(^)(t) 
t s u 

+ -eBub(6)(s) + -Tri(6)(t) + -Tri(6)(s) 

M U U 



St 



+ -^Box(^) (s, t) + - e(l - e) Box^^) (s, t) 



(3.6) 



u 

Here Bub(")(s), Tri(")(s) and Box(")(s,t) are the one-loop bubble, triangle and box 
scalar integrals, evaluated m. D = n — 2e dimensions. The remaining helicity configura- 
tions can all be obtained using parity and relabelings. In the CDR scheme one would 



need additional e-helicities 32 



The bubble and box integrals that appear in the above amplitudes are 

Bub(^)(s: 



e(l-2e) 



Bub(*^)(s) 
Tri(^)(s) 
Tri(6)(s) 



rr 



2e(l -2e)(3-2e) 

■5 (-)-'-. 



rr(-s) 



where 



rr 



2e(l -2e)(l -e) ' 



r(l + e)r2(l-e) 



(3.7) 



(3.8) 



r(l - 2e) 



(3.9) 



8 



with 



n 



n=l 



S,2 — — 
D 



C3 = 1.202057...^ 



(3.10) 



and we have kept the full dependence on e in the integrals. In the s-channel where 
s > the functions are given by the analytic continuation ( p.l2| ). 

The box integrals in various dimensions appearing in eq. (|3.6| ) are related via a 
on-shifting 

Box(^)(s,t) 



dimension-shifting formula ||33| valid to all orders in e, 

1 



Box(«)(s,t) 



2(-l + 2e)u 
1 

2(-3 + 2e)n 



stBox(''(s,t) - 2tTri(')(t) - 2sTri(')(s) 
stBox(*^)(s,t) -2tTri(^)(t) -2sTri(^)(s)) . (3.11) 



Since the D = Q — 2t scalar box integral is completely finite as e — 0, it is convenient to 
express the other box integrals in terms of it. This isolates all divergences to triangle 
and bubble integrals. In the u-channel (s < 0, t < 0), where the functions are manifestly 
real, the expansion of the six- dimensional box through 0{e^) is [llO| , p 



-l-e 



Box(^)(s,t) 



2(1 - 2e) 

+ 2t(hi,{-v) - VU^{-v) - — - 

\ 2 o o 

T 2 2 

+ tv^w^ - - ^—VW - 2C4) + (s ^ t) 



+ 0{e 



(3.12) 



where 



w 



r = ln - 



u 



u 



u 



W = \n - 



t 



u 



(3.13) 



In the s-channel (s > 0, t < 0) an analytic continuation of the box integral yields, 



Box(^)(s,t) 



rr s 



M(l-2e) 12' 
+ ef - Li3(-a;) +XLi2(-x) 



+ (3 + -YX^ - -Ti'^X 
3 ^ 2 2 



- e'[hu{-x/y) - hu{-y) + Li3(-y)X 

+ -U2{-x){X^ + TT^) + -iY^ + T,^f - -Y^X + -Y^X^ 

-r 2 A A T ' ' 6 4 
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+ Ix^r --iX^ + 7r2)2 + —XY + — TT^ 
3 8^ ^3 360 



+ m 



X + e(Li2(-a;) + FX - -X^ - ^ 



"6 

+ e\-U,{-x) - Li3(-i/) - \yX^ + + C. 



(3.14) 



where the variables appearing in the amphtudes are 



t 



u 



X 



X = In - 



t 



Y = ln 



(3.15) 



In the permutation sum appearing in eq. (p.2|), after appropriate relabehngs, both 
expansions of the box integrals ( |3.12| ) and ( ^.14] ) appear. 

Through C(e°) the amplitudes simply greatly and the M^^^^^^^^^ reduce to, 



M 
M 



(1) 

++- 

(1) 



M 



(1) 



1 + 0{e) , 

MiU. = Mil 

1 +2 I „,2 



M 



(1) 



1 + 0{e) 



M 



(1) 



M|'l_+(s,t,«) =Mi^j_+(s,n,t). 



2 s2 
2 m2 



(1) 













'-"in(') 









ln^(-^)-^lnf-^) 



in 



r + 



+ 0{e) 



(3.16) 



4. Finite parts of the two-loop gg 77 amplitudes 

A generic sample of two-loop Feynman diagrams for gg — ^ 77 is shown in Figure |^. 
We did not evaluate the diagrams directly. Instead we computed the unitarity cuts 
in various channels, working to all orders in the dimensional regularization parameter 
e = {4 — D)/2 |T^. Essentially we followed the approach first employed at two loops 
for the pure gluon four-point amplitude with all identical helicities and for X = 4 
super symmetric amplitudes ||^. These amplitudes were simple enough that a compact 
expression for the integrand could be given. The fermion loop contributions with all 
plus helicities are about as simple [Q. However, for the generic helicity configuration, 
the integrands become rather complicated. We therefore used the general integral 
reduction algorithms developed for the all-massless four-point topologies [|l^, ^ |19[, in 



order to reduce the loop integrals to a minimal basis of master integrals. To incorporate 



10 



polarization vectors of photons and gluons with definite hehcity requires some minor 
extensions of these techniques . 






Figure 2: Some of the two-loop diagrams for gg — > 77. The curly lines represent gluons 
while the wavy ones photons. 

We then expand the master integrals in a Laurent series in e, which begins at order 
Many of the master integral expansions quoted in refs. |T^, [T^, |T^, ^ are in 
terms of Nielsen functions [B6i , 



at 

/ _ln"-4 In^'fl-xt), 
! JO t 



(4.1) 



(n — l)!p! JO 

with n + p < A. We have found it useful to express the results instead in terms of a 



minimal set of polylogarithms IpTj , 

Li„(x) = £ ^ = /" ^ Li„_i(t) , 
Jo t 

Li2(x) = - / -ln(l-t) 
JO t 



(4.2) 
(4.3) 



with n = 2,3,4, using relations such as |38 



Si3{x) = - Li4(l -x)+ ln(l - x) Li3(l + - ln'(l - x)(Li2(x) - C2 



+ -ln^(l - a;) lna; + C4 



5*22(3;) = 114(3;) - 114(1 - x) + Li4 



-X 



1 — X 



ln(l - x)(Li3(x) - Cs 
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+ — ln^(l -x) 
for < X < 1. 



- ln=^(l - x) Inx + -C2 ln'(l - x) + (4 , 
2 



(4.4) 



The analytic properties of the non-planar double box integrals appearing in the 
amplitudes are somewhat intricate [1T0| , |T^, since there is no Euclidean region in any 
of the three kinematic channels for the 2^2 process. We quote our results in the 
physical s-channel (s > 0; t, u < 0) for the gg 77 kinematics ( ^.9| ). 

The dependence of the the finite remainder in eq. ( |2.14| ) on quark charges, the 
renormalization scale /i, and may be extracted as, 



_yV^{2)fin 

99^11 



26' 



ab 



-'A1A2A3A4 



2Nf 



6 



(ln(/iV. 



(1) 

Ai A2A3A4 



Ai A2A3A4 



N 



A1A2A3A4 



(4.5) 



where the spinor phases are defined in eq. (|3.5| ) and coupling constants have been 
extracted in eq. ( p.3|) . The /i-dependence in the first term in this expression is a 
consequence of renormalization group invariance. 

The two-loop QED corrections to gg — ^ 77 require the same set of two-loop di- 
agrams as the subleading-color QCD corrections. In the QED case, external fermion 
bubble insertions on the photon legs should be added. However, as mentioned in sec- 
tion ^ these diagrams are exactly cancelled by the ultraviolet counterterm in the usual 
on-shell renormalization of QED (just as in the light-by- light case |3^). The renormal- 
ized amplitudes are free of infrared and ultraviolet divergences and are given by 



A/f(2)QED 



.Nf 



'A1A2A3A4 -'^AiA2A3A4 5 



(4.6) 



where a factor of a(/i^) replaces one factor of as{jJ^) in the prefactor in front of AA^'^'^ 
in eq. ( p.3|) . The two- loop QCD and QED corrections to light-by-light scattering via 

(In ref. |^ particles 1 and 
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massless fermions are also proportional to -FAa2A3A4 
2 are taken to be incoming, so the helicity labels for legs 1 and 2 are reversed with 
respect to the labeling used here.) 

The explicit forms for the F^^^^^x^x^ appearing in eq. ( [4. 5] ) are 

_ 1 
~ 2' 
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(4.10) 



where the last amphtude does not possess symmetry under t ^ u. In contrast to 

the one-loop case, the — case is distinct from the +H h case, because at two 

loops external gluons can couple to a gluon internal to the diagram via the non-abelian 
coupling, while photons cannot. 

Similarly, the subleading color contributions in eq. ( [4.5|) are expressed in terms of 
the functions, 
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(4.15) 



The variables appearing in the amphtudes are defined in eq. ( |3.15|) . 



(4.16) 



5. Finite parts of the two-loop amplitudes 



The finite remainder for the (77 (77 amphtude is defined in eq. (|2.15| ). Its dependence 
on quark charges, the renormahzation scale /x, X and Nf is given by 
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(5.1) 



15 



In this case legs 1,3 are the gluons and legs 2,4 are the photons. The subleading-color 
functions -^aa2A3A4 same as for gg 77. 

The leading-color functions G\_^)^^)^^^^ are new, and are given by 
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where again the variables appearing in the amphtudes are defined in eq. ( |3.15D . 

In comparison with the one-loop amplitudes for (77 — > (77, the two-loop amplitudes 
are more singular in the forward limit m — > by one power of s/u, due to the exchange 
of a pair of gluons in the w-channel. The power-enhanced terms are only found in the 
leading color functions G^, and only in the helicity configurations for which the helicity 
of the incoming gluon does not fiip as it scatters (in the conventional helicity labeling, 
not the all-outgoing one used here). The helicity of the photon may or may not fiip in 
such terms. Explicitly, the power-enhanced terms are 



0, 
2m- 



u 



9 u 
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~ — (vr^ + 3)i7r— . 
9 u 



(5. 



These terms represent the beginning of the Reggeization of the gj —>■ gj amphtude in 
perturbation theory; the large logs of ln(— s/m) that will arise at subsequent orders in 
as could be resummed using BFKL techniques EO . 



6. Checks on results 

We performed a number of consistency checks on the amplitudes to ensure their relia- 
bility: 

1. As a check of gauge invariance, we verified that the amplitudes vanish when a 
gluon or photon polarization vector is replaced with a longitudinal one. 

2. The agreement of our explicitly computed infrared divergences with the expected 
form provides a stringent check on the amplitudes. Since the integrals generally 
contain both divergent and finite terms, this also provides an indirect verification 
that the leading color finite remainders have been assembled correctly. 



3. Using Supersymmetry Ward Identities we evaluated the identical helicity 
case by relating it to the already known identical helicity four-gluon ampli- 



tudes |T0[- Since the integration was done by a completely different technique, 
the agreement between the two independent ways of evaluating this amplitude 
provides an additional stringent check of the programs and integration methods 
used to obtain the general helicity cases. 

We compared the results for gg gg obtained using our methods and computer 
programs |^ to those of ref. ||13|. The interference of the two-loop gg —>■ gg 
helicity amplitudes with the tree amplitudes, after summing over all external 
helicities and colors and accounting for the different scheme used (HV vs. CDR), 
is in complete agreement with the calculation using conventional dimensional 
regularization ITH . 



7. Conclusions 

In this paper we presented the two-loop matrix elements for two gluons to scatter into 
two photons. Due to the large gluon-gluon luminosity in the x range relevant for the 
LHC Higgs search, these two-loop contributions should be competitive in size with the 
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existing next-to-leading order corrections to the qq — > 77 subprocess 0. The matrix el- 
ements presented here, together with the associated real emission contributions, namely 



the one-loop amplitudes for gg — > 77(7 p4|, ^ , will feed into improved estimates of the 



QCD background to Higgs production at the LHC, when the Higgs decays into two 
photons 0. 
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